Charmonium in finite temperature lattice QCD 
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We study hadron properties near the deconfining transition in the finite temperature lattice QCD. This paper 
focus on the heavy quarkonium states, such as J/i/; meson. We compare the meson correlators above and below 
Tc and discuss the possibility of the cc bound state by observing the wave function. 



1. Introduction 

It is widely believed that the Quantum Choro- 
modynamics (QCD) exhibits a phase transition 
at some temperature Tc, and quarks and gluons 
confined in the low temperature phase are liber- 
ated to form "quark gluon plasma". At the be- 
ginning of 2000, CERN reported that the QGP 
state had been created in the heavy ion collision 
experiment ||l[. In this experiment, the J/tp sup- 
pression is regarded as a key signal of QGP 
formation. More reliable inspection will be per- 
formed in RHIC project at BNL. 

On the theoretical side, in spite of various ap- 
proaches, we are still far from definite under- 
standing of charmonium properties near the tran- 
sition and the fate of the hadronic states in the 
plasma phase. We investigate this problem using 
lattice QCD, which enables us to incorporate the 
nonperturbative effect of QCD. 

2. Our approach 

The difficulties in studying the heavy quarko- 
nium state at finite temperature can be catego- 
rized into two classes: (i) difficulties to obtain de- 
tailed information on the heavy quarkonium cor- 
relators at finite temperature, and (ii) analyses of 
the obtained correlators. Although it is difficult 
to overcome the latter, the former is possible to 
be solved by employing the techniques developed 
in recent years. 

At high temperature, sufficient information on 
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the temporal correlator requires high resolution in 
the temporal direction. To achieve large tempo- 
ral cutoff with limited computational resources, 
we employ the anisotropic lattice, on which the 
temporal lattice spacing is finer than the spatial 
one 1^. Calculations with a heavy quark need 
further implementation. The standard lattice ac- 
tion contains large systematic uncertainty for the 
quark mass larger than or comparable to the lat- 
tice cutoff. This problem was reexamined with 
the light of effective theoretical interpretation of 
the Wilson (and clover) quark action by FNAL 
group |4|. Their approach is well-matched with 
to the anisotropic lattice, since it naturally in- 
troduces the anisotropy in the spatial and the 
temporal hopping parameters. We employ this 
approach to describe the heavy quark on present 
lattices, which have not sufficiently large spatial 
cutoff compared with the charm quark mass. 

The latter problems, (ii), are much more diffi- 
cult to be overcome. This is because one is in- 
evitably enforced to extract the spectral proper- 
ties from the data at the short temporal distance, 
where high frequency component of the dynam- 
ics is significant. At present, there is no well- 
established way of lattice simulations to attack 
the spectroscopy at T > 0. For this reason, lat- 
tice studies of hadrons at finite temperature have 
been argued on the spatial correlation (screening 
mass) [|. Recently, QCD-TARO CoUaboration 
(including three of us) analyzed the space-time 
structure of the mesonic correlators in the tempo- 
ral direction to examine the temperature effect on 
them They caught a sign of the bound state 
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even in the deconfined phase, as well as an indi- 
cation of the chiral symmetry restoration. Their 
strategy seems one of the best approaches also to 
the present subject, and is expected to give us 
novel information on the fate of charminium near 
the phase transition. Therefore we proceed the 
analysis with the following steps. 

• At T = we construct good meson operators 
which has large overlap with the state of inter- 
est. 

• Temperature effects are observed in the meson 
correlator between the operators constructed in 
the previous step. 

• We define the "wave function" in the Coulomb 
gauge at finite temperature, and observe its t- 
dependence below and above Tc to discuss the 
possibility of the bound state at T > Tc. 

3. Fermilab action on anisotropic lattice 

The quark action we use is in the following 
form. 

Sf =Y.q{x)K{x,y)q{y), (1) 

-Kr{(l - 74)r+4 + (1 + 74)^-4} 

-^■yY.iii'^ - li)T+i + {r + lt)T^i} 

-{KaCEgo-iiFiiix) + rK^CB-gcrijFij{x))Sx,yA'^) 

where T±^ = U±fj_{x)Sx±p,.y The bare anisotropy 
of the fermion field appears in the ratio of 
the spatial and the temporal hopping parameters, 
Kt = "fpKa- The Wilson parameter in this form is 
r = 1/^, where ^ is the renormalized anisotropy: 
^ = aa/ur- The hopping parameters are related 
to the bare quark mass as Ko- = l/2{mQ+jF+3r). 
We also define n with 1/k — 1 / k„ — 2{jp + 3r ~ A) 
as the parameter which plays the same role as k 
on the isotropic lattice. As the field strength F^^,, 
we use the standard clover-leaf definition. At the 
tree level, the coefficients of the clover terms, ce 
and CB, are unity. 

This action is obtained by generalizing the form 
in Q to the anisotropic lattice, and by setting 
r — 1/^. This choice of the Wilson term is 
adopted so that the action restores the axial inter- 
change symmetry in the physical unit as niQ — > 0. 
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Table 1 

The quark parameters and the meson masses in 
the physical unit. The error of the meson masses 
include the statistical error of the scale a^^. 



Therefore we can use the same form in the light 
quark region as the 0(a) improved Wilson quark 
action . This is the reason that our form of the 
action is different from the form used in |^,^ . 

On the anisotropic lattice, we need additional 
tuning of the parameters called as "calibration" , 
to ensure the renormalized anisotropy is same for 
the gauge and the fermion fields [|j . We define the 
fermionic anisotropy S^p using the dispersion re- 
lation of the pseudoscalar and the vector mesons 
and tune the value of 7^? so that ^p coincides with 
the gauge field anisotropy ^. 

We incorporate the mean-field improvement 
0: U, -> Ujuoa (i =1,2,3) and C/4 -> U^/uor 
with the mean-field values of the spatial and the 
temporal link variables, uq^ and uqt defined in 
the Landau gauge. 

4. Result at zero temperature 

Numerical calculation is done on the lattice of 
sizes 12^ X Nt where Nt = 72, 20, 16 and 12, with 
the standard Wilson gauge action at the coupling 
(3 = 5.68 and the bare anisotropy 70 = 4.0, in 
the quenched approximation ^ . On Nt = 72 lat- 
tice, the renormalized anisotropy and the spatial 
lattice spacing from the string tension are deter- 
mined as C = 5.3(1) and = 0.85(3) GeV. 
Nt = 20, 16 and 12 lattices correspond to the 
temperatures T ~ 0.93 T^, 1.15 and 1.5 T^, 
respectively. We calculate the meson correlators 
at r = with the three sets of the parameters 
(k, 7f) listed in Table ^ 

Each value of jp is determined by the calibra- 
tion procedure mentioned above. Observing the 
resultant meson masses summarized in Table |l|, 
(k, jp) = (0.1013, 3.703) roughly corresponds to 
the charm quark. In the successive analysis, we 
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Figure 1. Effective mass plot for the diagonalized 
correlators in the vector channel. 
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Figure 2. Effective mass plot at finite tempera- 
tures in the vector channel. 



use this set of parameters. 

To optimize the meson correlator at T = 0, we 
employ the variational analysis. This is imple- 
mented by diagonalizing the correlator matrix 



(3) 



The operator Oi (x, t) is defined with the smearing 
function ip in the Coulomb gauge as 



(4) 



where the 4x4 matrix F specifies the quantum 
number of the meson. As the smearing function 
ipi, we use the eigenfunction obtained by solving 
the Schrodinger equation for the S state, 
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(5) 



where V{r) is the static quark potential obtained 
on the same lattice, and m/j = 1.5/2 [GeV] is 
used. 

Figure |l] shows the effective masses of the di- 
agonalized correlators with two largest eigenval- 
ues (smallest masses). Both correlators show 
plateaus beyond t ~ 10, though the statistical 
fluctuation is large for the correlator correspond- 
ing to the first excited state. We fit these data 
aX t — 16-36 and find the meson masses in 25 
states as mps[2S) — 0.7677(77) and TOy(2S) — 
0.7794(88) in the lattice unit. Although these 



give 15-25 splitting smaller than the experimen- 
tal value, considering the large systematic uncer- 
tainty on this coarse lattice, they are not inconsis- 
tent results. We regard the operator corresponds 
to the ground state in this analysis as the opti- 
mized meson operator at T = 0. 

5. Result at finite temperature 

We start the argument at T > with the tem- 
perature dependence of the correlator between 
the optimized operators at T = difined in the 
previous section. Figure |^ shows the effective 
mass plot at finite temperatures in the vector 
channel. Below Tc {Nt = 20), the plateaus of 
the effective mass seem to appear. Although the 
effective mass of the ground state is slightly larger 
than that at T = 0, it is difficult to identify the 
plateau precisely and determine the mass quan- 
titatively with present statistics. More detailed 
analysis with higher statistics may open a stage 
to discuss the potential mass shift of charmonium 
near to the 

Above Tc {Nt = 16 and 12), the correlators do 
not show any clear plateau in the effective mass 
plot. They decrease more rapidly than expected 
from the behavior at Nt — 20 in both of the cor- 
relators correspond to the ground and the first 
excited states at T = 0. The Ps channel shows 
similar feature, although the decrease of the ef- 
fective mass is slightly milder. These behavior 
at least signal significant change in the nature of 
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Figure 3. The t-dependence of the wave func- 
tion normahzed at the spatial origin in the vector 
channel at Nt =20 and 12. For comparison, the 
same quantity composed of the free quark propa- 
gators is also shown. 

the correlators when the system crosses Tc- In ||], 
the (effective) masses increase as T in the pseu- 
doscalar and the vector channels. The observed 
behavior in present work, however, shows quali- 
tatively different nature of the correlators. 

To discuss the spatial correlation between 
quark and antiquark at T > 0, we observe the 
"wave function" normalized at the spatial origin, 

0(r,t) =ii;r(r,<)/wr(0,i), (6) 
^T{r,t) = Y.M^+r,t)Vq{x,t)0\Q)). (7) 

If there is no bound state, like as with the 
free quark propagators, this wave function (/)(r, t) 
broadens with t for any source smearing function. 
In contrast, the existence of the stable shape of 
(j){r, t) give us a hint on the existence of the bound 
state. 

Figure || shows the i-dependence of the wave 
function normalized at the spatial origin in the 
vector channel at Nt =20 and 12. The source 
function of these correlators are taken to be wider 
than the optimized one, to show clearly that 
4>{r,t) keeps narrow shape. At the both Nt, the 
wave function gradually decrease as t, in con- 
trast to the case with the free quark poropagators 
shown together. This indicates that the quark 
and antiquark tend to stay together even in the 
deconfined phase, up to 1.5 T^. 



These two observations are interesting, and at 
the same time puzzling. Obviously, the phase 
transition causes significant change in the char- 
monium correlators, especially in the vector chan- 
nel. On the other hand, the result of the wave 
function suggests the persistence of the bound 
states. One natural picture is that the mesonic 
spectral function still have a peak with rather 
large width above Tc. This explains both of our 
observations. Direct extraction of the spectral 
function from the lattice data jl2| may give us fur- 
ther information. For the definite understanding 
of the fate of chamonium at the phase transition, 
we need more systematic studies and development 
of procedures. 

The simulation has been done on Intel Paragon 
XP/S and NEC HSP at INS AM, Hiroshima Uni- 
versity. This work is supported by the Grant-in- 
Aide for Scientific Research by Monbusho, Japan 
(No.l0640272,No.ll440080) 
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